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Beamforming on the MISO interference channel 
with multi-user decoding capability 
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^ ' This paper considers the multiple-input-smgle-output interference channel (MISO-IC) with in- 

terference decoding capability (IDC), so that the interference signal can be decoded and subtracted 
from the received signal. On the MISO-IC with single user decoding, transmit beamforming vectors 
are classically designed to reach a compromise between mitigating the generated interference (zero 



ry^ ' forcing of the interference) or maximizing the energy at the desired user. The particularly intriguing 

problem arising in the multi-antenna IC with IDC is that transmitters may now have the incentive 



to amplify the interference generated at the non-intended receivers, in the hope that Rxs have a 
better chance of decoding the interference and removing it. This notion completely changes the 
previous paradigm of balancing between maximizing the desired energy and reducing the generated 



^3 I interference, thus opening up a new dimension for the beamforming design strategy. 



Our contributions proceed by proving that the optimal rank of the transmit precoders, optimal 
in the sense of Pareto optimality and therefore sum rate optimality, is rank one. Then, we inves- 
tigate suitable transmit beamforming strategies for different decoding structures and characterize 
the Pareto boundary. As an application of this characterization, we obtain a candidate set of the 
maximum sum, rate point which at least contains the set of sum rate optimal beamforming vectors. 
We derive the Maximum-Ratio- Transmission (MRT) optimality conditions. Inspired by the MRT 
optimality conditions, we propose a simple algorithm that achieves maximum sum rate in certain 
scenarios and suboptimal, in other scenarios comparing to the maximum sum rate, p 

This work has been performed in the framework of the European research project SAPHYRE, which is partly 
funded by the European Union under its FP7 ICT Objective 1.1 - The Network of the Future. This work is also 
supported in part by the Deutsche Forschungsgemeinschaft (DFG) under grant Jo 801/4-1. 

^This work was partially presented in [1]. 
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I. Introduction 

Despite efforts since pioneering work such as [2], [3], the capacity region of interference channel 
is still an open problem. Numerous work have attempted to compute achievable rate regions and 
outer bounds on the Single-Input-Single-Output Interference Channel (SISO-IC). It is proved in [2] 
that the capacity of any two-user SISO-IC is the same as in a corresponding IC in standard form: 
direct gain as unity and interference gain as a real positive scalar. Even since, the capacity region 
of two-user SISO-IC has been studied extensively (see e.g. [4]-[8] and the references therein) : 

• In the weak interference regime where the cross interference gain is much weaker than the direct 
channel gain, the sum rate capacity is shown to be achievable by treating interference as thermal 
noise at the receiver which requires no feedback communication between Rx j and Tx i [4] , [5] . 

• At the other extreme, in the strong and very strong interference regime, both users should decode 
the interference signal while treating the desired signal as noise. The decoded interference signal 
is then subtracted from the received signal allowing the desired signal to be decoded without 
any interference. [2], [9], [10]. 

• In the mixed interference regime, where one cross interference gain is stronger than direct channel 
gain and the other link is weaker, the sum rate capacity is shown to be attained by one user 
decoding interference and the other user treating interference as noise [5], [8]. 

• The deterministic channel approach offers a good approximation of the sum capacity of inter- 
ference channel. In the deterministic channel approach, the input-output relationship of the 
channel is modeled as a bit-shifting operation [7], [11], [12]. In [7], the two-user SISO-IC sum 
capacity is approximated to within one bit using the deterministic channel approach. 

To extend the above results, the conditions in which treating interference as noise achieving 
capacity on the vector Gaussian interference channel is studied in [13]. The capacity region of a 
specific class of MIMO interference channels is characterized in [14]. Against intuition, the optimal 
rank of input covariance matrices remains inconclusive, unlike in single user detection (SUD) case 
where single mode beamforming attains capacity [4]. The authors in [14] showed that the optimal 
input covariance matrix attaining capacity of MISO-IC has rank less than the number of users in 
the IC. 

The frontier of the achievable rate region, also known as the Pareto boundary, holds importance 
to the understanding of IC. Any rate points on the Pareto boundary are operating points such that 
one user cannot increase its rate without decreasing other users rates. Assuming perfect CSIT, the 



Pareto boundary of SISO-IC and MISO-IC with SUD are characterized in [15], [16] respectively. 
In [17], the authors extended the results to partial CSIT. In this paper, we assume simple single 
user encoding transmitters and interference decoding capability at receivers, which yield a simpler 
scheme comparing to the Han-Kobayashi scheme [9] . We study the effects of transmit beamforming 
on the achievable rate region and to characterize the Pareto boundary. We limit ourselves to the 
two transmitter-receiver (Tx-Rx) pairs interference channel with IDC. We assume each receiver 
can choose to fully decode interference (D) or treat interference as noise (N). No rate splitting- 
based averaging is considered between these two modes. The transmit beamforming vectors are 
optimized to achieve an operating point as close to the Pareto boundary as possible. In IC with 
SUD, interference mitigation may seem to be a reasonable strategy. Yet, with the IDC which we 
address in this paper, it is possible for the user to manipulate it's beamforming vector such that the 
generated interference is amplified for easier interference removal and yield a better operating point. 
The fundamental question becomes: when should the beamforming vectors be designed to amplify 
interference to improve performance and when to mitigate interference? 
The main contributions of this paper are: 

• In Section IIIII we describe an achievable rate region of the MISO-IC with IDC, with the 
assumption of linear precoding, taking into account of receivers choice of actions, D or N. 

• We study and characterize its Pareto boundary in Section lVl-A| in terms of beamforming vectors 
design and power allocation. We characterize the set Q of tuples of beamforming vectors and 
power allocation which attain the Pareto boundary. 

• As a special case, in Section lVlH we characterize the set of beamforming vectors which attain the 
maximum sum rate point in the form of a candidate set Q. As the maximum sum rate problem 
is non-convex, conventional solutions rely on different searching techniques. Note that Q G fl. 
The cardinality of Q is much smaller than the cardinality of Q which provides a significant 
reduction of searching complexity. Further, we prove that with IDC full power must be used 
at each Tx to attain the maximum sum rate point. This result is interesting as non-full power 
should be employed in some Txs to achieve the maximum sum rate point in the SISO-IC-SUD 
[15], [18]. 

• In Section IVIIH we investigate the conditions of channel parameters for which simple strategies 
are sum rate optimal. In particular, we study the matched filter (MF) with respect to the 
desired channel and the MF with respect to the interference channel, which are termed as the 



maximum-ratio-transmission (MRT) schemes. 

• Inspired by the MRT optimahty conditions, we propose a suboptimal but very low complexity 
beamforming design in Section IIX-DI The suboptimal algorithm shows a promising tradeoff 
between complexity and performance, as illustrated by simulation results. 

• In Section HXl we provide simulations and discussions which illustrate cases where interference 
decoding is most beneficial to sum rate performances. 

A. Notations 

The lower case bold face letter represents a vector. The conjugate transpose is denoted by (.)^. 
M represents the set of real numbers. The projection matrix on vector x is IIx = xx^/||x|p and the 
orthogonal projection matrix is I — IIx where I is the identity matrix. Denote a boolean statement by 
Bi. The complement of the statement Bi is Bi. The OR operation is denoted as U; AND operation 
as n. i'{A) returns the dominant eigenvector of matrix A. tr{A) is the trace of matrix A. The 
matrix A is positive semi-definite if A ^ 0. The symbol <^ represents the if-and-only-if relationship 
between two statements. Re{z) and Im{z) give the real and imaginary part of complex number z. 
The function arg(z) gives the phase of the complex number z. The operator x is the Cartesian 
product operator between two sets. 

II. Channel model 

We assume a system of two transmitter-receiver (Tx-Rx) pairs in which each Tx has A^ transmit 
antennas and each Rx has only one receive antenna. This results in a two-user Multiple-Input-Singie- 
Output Interference Channel (MISO-IC), which is illustrated in Fig. [T]as an example with A^ = 3. 
We assume that the Txs are using commonly known codebooks and therefore the Rx, if the channel 
qualities allow, can decode the interference and subtract it from the received signal. Also, we assume 
that the interference is successfully decoded if the rate of the interference signal is smaller than the 
Shannon capacity of the interference channel. 

In the MISO-IC-SUD, it has been shown that the optimal transmit precoders are rank 1 and 
therefore beamforming attains the Pareto boundary. However, whether his conclusion holds in the 
MISO-IC-IDC is not known yet. We answer this question in the following by starting with a general 
transmit covariance matrix. Denote the transmit covariance matrix of Tx i by S^ and the channel 
from Tx i to Rx i, where i £ {1,2} ,i 7^ i, h^^ E C^^^. Note that the channel gains are proper 
i.i.d complex Gaussian coefficients with zero mean and unit variance. The received signal at Rx i is 




Fig. 1: The 2 users MISO-IC where Txs are equipped with 3 antennas. 
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The noise rii is a complex Gaussian random variable with zero mean and unit variance. The symbol 
Xi is the transmit symbol at Tx i with unit power. Denote the set of the transmit covariance matrices 
that satisfy the power constraint ir(Si) < Pmax to be 

^NxN 



5 = {S G C^^^ : S ^ 0, tr(S) < P^a.) , i = {1, 2}. 



(2) 



in. Achievable Rate Region 

We propose the following four decoding structures corresponding to the Rxs. actions: (N,N), (N,D), 
(D,N) and (D,D) [19], with "Nf^tands for treating interference as noise and "Df^tands for decoding and 
removing interference. Thus, (D,N) means Rx 1 decodes and removes interference and Rx 2 treats 
interference as noise. In [19] , these four decoding structures are proposed and its corresponding rate 
points are shown to be achievable in the SISO-IC. We extend the concept to the MISO-IC and define 



the following important quantities: 

Ci(Si)^log2(l + hfiSihn), 

C2(S2)^log2(l + hiS2h22), 

' hfiSihn 



i^l(Si,S2)^log2ll 
i^2(Sl,S2)^log2fl 
T2(Si,S2)^log2(l 



hf2S2hi2 + i; ' 

_h|S2h2^\ (3) 

hfiSih2l + lj' 

hi2S2hi2 \ 

hfiSihn + 1/' 



^■(^'•«^)*'-^t-ii£TT) 



C\ and C2 are the single user rates, the largest rate user 1 and 2 can achieve without the influence 
of interference. Di and D2 are the rates corresponding to decoding the desired signal while treating 
interference as thermal noise and Ti and T2 are the rate corresponding to decoding the interference 
while treating the desired signals as noise. 

Consequently, if both receivers decode interference, user i must transmit at a rate that ensures 
interference decoding at Rx i, thus we have the following: 

i?i<min{Ci(Si),Ti(Si,S2)} 

(4) 
i?2<min{C2(S2),T2(Si,S2)}. 

Denote the rate region with interference decoding at both receivers by the Decode-Decode (DD) 

region: 

n''''= U |(i?i,i?2)< (min{Ci(Si),Ti(Si,S2)},min{C2(S2),r2(Si,S2)})|. (5) 

(Si,S2)g5xS *- ^ ^ ^ 

Remark 1: For each selected pair of transmit beamformers, a corresponding rate region which 

satisfies the inequalities (JH) is obtained. The achievable rate region TZ'^'^ is defined as the union of 

all regions achieved by all possible transmit beamformers. 

On the other hand, if both Rxs choose to treat interference as noise, we obtain the NN region, 

7^™= U |(i?l,i?2)<(l?l(Sl,S2),D2(Sl,S2))|. (6) 

If Rx 1 decodes interference but Rx 2 treats interference as noise, Tx 2 must transmit at a rate 
that ensures interference decoding at Rx 1. Thus, the DN region is obtained as, 

n''^= IJ |(/?i,i?2)< (Ci(Si),min{i^2(Si,S2),T2(Si,S2)})|. (7) 



Remark 2: 7^'^"(Si,S2) is the rate region that the inequahties in ([T]) are satisfied for specific 
transmit covariance matrices (Si, 82)- It can be an empty region if the inequahties cannot be satisfied 
at the same time. This corresponds to the situation where the data rate of Tx 2 is too high for Rx 
1 to decode. 

Similarly, we have for the ND region, 

n^'^= U |(i?i,i?2)< ('min{A(Si,S2),Ti(Si,S2)},C2(Si))|. (8) 

Finally, one achievable rate region of the MISO-IC with interference decoding capability is therefore 
the union of the above regions: 

^ ^ ^nn ^ -j^dd ^ -j^dn ^ ^nd _ ^g^ 

We now turn our attention to the Pareto boundary of the rate region. To find the boundary 
achieving solutions, we proceed by identifying a set of smaller dimension than S x S which is 
guaranteed to contain the Pareto optimal solutions. We refer to such set as candidate set. The main 
practical value of a candidate set is to offer a subtantial reduction of complexity compared with the 
exhaustive search over the full set S x S. 

Definition 1: Denote the set of points on the Pareto boundary of TZ by B{TZ). If the rate pair 
(ri,r2) € 7^ is on the boundary, (ri,r2) E B{TZ), then there does not exist a rate pair (?'i,T2) G TZ 
such that (^'^,^2) > (ri,r2), with one strict inequality. Using TZ in ^, 

B{TZ) c B(7^"") u BiTZ'^'^) u B(7^'^") u B(7^"'^) (lo) 

Definition 2: The transmit covariance matrices Si, S2 are Pareto optimal in the rate region TZ if 



i?i(Si,S2),i?2(Si,S2)j G5(7^). (11) 

Definition 3: The candidate set Q.^^ of B{TZ^y), a;,y G {n, d}, is a set of transmit covariance 
matrices that contains the transmit covariance matrices that attain the Pareto boundary of TZ^^ . If 
(Si,S2) are Pareto optimal, then (Si,S2) G ^'^^ . Similarly, the candidate set of B{TZ) is 17 which 
contains all pairs of (Si, S2) which are Pareto optimal in the region TZ. 

IV. The Pareto optimal transmit covariance matrices 

In this section, we study the transmit covariance matrices that attain the Pareto boundary and 
prove that they are rank one. 



Theorem 1: The Pareto boundaries of the NN region, the DN region and the DD region are 
attained by rank 1 matrices. Consequently, the Pareto boundary of MISO-IC-IDC, defined in (|10|) . 
is attained by rank one transmit covariance matrices, or transmit beamforming. 

Proof: Here, we provide a sketch of the proof. For details, please refer to Appendix IXI-AI We 
first show that the boundaries of rate region T^'"^ and TZ'^'^ are attained by rank one matrices. By 
exchanging the roles of the transmitters, we obtain that B{TZ'^"') is attained by rank one matrices. 
From [4], [20], it is shown that the boundary in the NN region is attained by rank one matrices. Hence, 
the boundaries of all decoding structures are attained by rank one transmit covariance matrices. Since 
the Pareto boundary of the proposed achievable rate region in the MISO-IC-IDC, defined in (|1U|) . 
is a subset of the union of the above boundaries, we conclude that this Pareto boundary is attained 
by rank one transmit covariance matrices, or transmit beamforming. ■ 

From Theorem [TJ we have established that the Pareto boundary is attained by transmit beam- 
forming vectors. To facilitate the following discussions, we define the transmit beamforming vectors 
Wi and transmit power Pi, for i = 1,2, 

S, = Wiwf P, (12) 

with ||wip = 1. As an abuse of notation, we write S as the set of all possible beamforming vectors, 

Wie5, 5 = {wgC^^^ : ||w|| = l}. (13) 

Consequently, we redefine the candidate sets in terms of transmit power allocations and beamforming 
vectors. The candidate set 17^^ oi B(TZ^y),x,y £ {n,d} contains the Pareto optimal beamforming 
vectors and transmit power allocations. 

n^y D |(wi,W2,Pi,P2) : (Pi(wi,W2,Pi,P2),i?2(wi,W2,Pi,P2)) eS(7^^^)| (14) 

and the candidate set of B{TZ) is 

C!d|(wi,W2,Pi,P2): (Pi(wi, W2, Pi, P2), P2(wi, W2, Pi, P2)) G B(7^)|. (15) 

In the following sections, we study the Pareto boundary in terms of power allocation and transmit 
beamforming vectors in different decoding structures namely R"** and R'^'^. R"" is the case of MISO- 
IC-SUD and is well studied in [16]. R'*" is symmetric to R"'^ and is therefore omitted here. Then 
as a special case, we discuss the characterization of the maximum sum rate point in each decoding 
structures. 



|hj(wj2/^ 




l,0< R<P„ 



luH^, |2p. 



Fig. 2: The received power region for Tx i. 



V. The received power region 

The received power region was first proposed in [20] as a powerful tool to illustrate the dependancy 
between the received power tuple and the Pareto boundary on the i^-user MISO-IC-SUD. The tuple 
at one receiver includes the received power from the desired signal and the received power from 
the interference signal(s). In a two-user MISO-IC, we can illustrate the received power region as a 
two-dimensional plot, as shown in Fig. [21 Mathematically, the received power region of user i for the 
two-user MISO-IC is defined as: 



$, = {(|hf W,|2P,, |hjfw,pP,) : W, G 5,0 < F, < Prr^aa] 



(16) 



where the desired channel power of user i is Ih^w^pPi and the interference channel power of user i 
is |hj^w,|2Pi. 

The importance and relevance of the received power region can be summarized in the following. 

• The boundary of received power region and the Pareto boundary: The received power region is a 
convex and compact region with respect to the received power values. The Pareto boundary of 
the MISO-IC-SUD with linear pre-coding, the NN region here, is shown to be attained by the 
received power values on the boundary of the received power region [16]. 

• Monotonicity of rates: The rate metrics defined in ([3]) are either monotonically increasing or 
decreasing with the channel powers. The optimization of such rates can be simplified by first 



computing the optimized channel powers and then the corresponding beamforming vectors that 
achieve such channel powers. 
Now we define the received power region achieved by beamforming vectors inside the Pareto 
boundary candidate set, Q, 

<P{n) = |(^|hfiWipPi,|hfiWi|2Pi,|hf2W2pP2,|hf2W2pP2) : (wi , W2, Pi, P2) € ^y (17) 

Immediately, we have the following relations: the Pareto boundary candidate set achieves a received 
power region which is a subset of the Cartesian product of the received power regions for Rx 1 and 
2, <I>i and $2, 

$(fi) C <^i X $2- (18) 

The beamforming vectors and power allocations in $1 x <I>2 contribute to the whole achievable rate 
region whereas the tuples in $(f^) only attain the Pareto boundary. This means that if we know 
$(17), we can achieve the Pareto boundary without searching over all beamforming vectors in the 
remaining space in $1 x $2- This reduces the search space from <I'i x <I'2 significantly. 

In the following sections, we compute 17"", il'^", f]"'', il'^'* which are candidate sets of the Pareto 
boundary of the corresponding regions: NN, DN, ND and DD. We define a candidate set of the 
overall Pareto boundary, Q, as the union of the candidate sets mentioned above: 

n= \J n^y. (19) 

x,y€{n,d} 

VI. Pareto boundary characterization 

A. Pareto boundary characterization in the ND region 

With decoding structure R"**, Rx 1 treats interference as noise and Rx 2 decodes and subtracts 
the interference signal from the received signal before decoding the desired signal. 
Theorem 2: The Pareto boundary B{TZ^'^) is attained by candidate set Q^'^ 

n^'i = { Wl , W2 , Pi = P^a.. < P2 < P^a.} (20) 

where Wi , W2 defined in (|2ip , are sets of beamforming vectors composed of linear combinations of 
two channel vectors; to attain the Pareto boundary, Tx 1 transmits with full power Pmax whereas 
Tx 2 transmits with less than full power P2 < Pmax- 

Proof: See Appendix IXI-BI ■ 
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W. = |w, : w, = VaI^^ + yi^A-i^ ; < A, < l| , ^,J = l,2,z ^ j. (21) 



In the ND region, Rx 1 treats interference as noise and Rx 2 decodes interference. As described 
by Thm. El the Pareto optimal transmit power for Tx 1 is to transmit at full power Pmax and less 
than full power for Tx 2. The interpretation is that Tx I's transmit power does not affect the rate 
performance of Rx 2 as the interference from Tx 1 is decoded and removed. On the other hand, Tx 
2 is not advised to transmit at full power because its increase of power will increase the interference 
power at Rx 1 and hence reduce the achievable rate of Rx 1. 

The Pareto optimal transmit beamforming vectors are parameterized in the sets Wi , yV2 as positive 
linear combinations of two orthogonal vectors. These two vectors are the desired channel projection 
onto the span and the null space of the interference channel. As shown in Fig. [3l the vectors in Wi 
and W2 are represented by blue regions. The blue regions cover from the point of zero interference 
power (Point A in Fig. [3]) to the point of maximum desired channel power (Point B) and the point 
of maximum interference power (Point C). A/[oving from point A to B and C on the Pareto boundary 
in Fig. El the interference power increases monotonically. On the left figure of Figure El we show 
the received power region of Tx 1, the channel powers between A and B correspond to a strong 
desired channel power of Tx 1 and a relatively small interference channel power from Tx 1 to Rx 2. 
These points may attain the Pareto boundary if the desired channel power of Tx 2 is weak. On the 
other hand, if Tx 2's desire channel power is large, the interference power from Tx 1 to Rx 2 must 
be increased to increase to interference rate Ti which limits rate rate of Ri as Ri = min(Ti,L'i). 
This is a novel concept comparing to the conventional single user decoding interference channel, the 
increase in interference power here is beneficial as it facilitates interference decoding and removal. 

In the next section, we investigate the Pareto boundary attaining beamforming vectors in the DD 
region. 

B. The Pareto boundary characterization in the DD region 

With decoding structure R*^**, both Rx's decode interference. The Pareto boundary attaining 
solutions are: 

Theorem 3: The Pareto boundary B{V/'^) is attained by candidate set 

n'''' = {Wi G Vl, W2 G V2, < Fi, F2 < Pma.} , (22) 

11 
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f]"^ = {Pi = P^,,, wi G Wi, W2 G W2, < P2 < F™,,} 



|h22W2| P2 



Fig. 3: The graphical ihustration of the candidate set 17"*^ as the shaded area which is a subset of 
the received power regions ^i. 



with Pareto boundary attaining beamforming vectors composed of two orthogonal channel vectors, 
specifically for i = 1,2, 



Vi 



w G 5 : w 



A,;^^ + VT^Xl-^^^O < \ < 1 



in-'h--i 



|n,ih,,| 



(23) 



and both Txs transmit at less than full power. 

Proof: See Appendix IXI-CI ■ 

Remark 3: Note that Wi in Thm.[2]and Vi defined here are different candidate sets. In particular, 
Wi is a set of vectors that are the positive linear combinations of Hjiha and Hj^hn whereas Vi is a 
set of vectors that are the positive linear combinations of Hahji and n^hj^. This difference is shown 
graphically in Fig. [3] and |31 

In the DD region, both Tx 1 and 2 decode interference and their choice of actions are symmetric. As 
described by Thm. [3l the Pareto optimal transmit power for Tx 1 and 2 are to transmit at less than 
full power. The Pareto optimal transmit beamforming vectors are parameterized in the sets Vi, V2 
as positive linear combinations of two orthogonal vectors. These two vectors are the interference 
channel projection onto the span and the null space of the desired channel. As shown in Fig. [H the 
vectors in Vi and V2 are represented by blue regions. The blue regions cover from the point where 
the point of maximum desired channel power (Point B) to the point of maximum interference power 
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Fig. 4: The graphical illustration of candidate set Q.'^'^ as the shaded area which is a subset of the 
received power regions ^i. 



(Point C) and the point of zero desired channel power (Point D). 

Notice that the points where the interference channel powers are zero is not Pareto optimal. It 
is because minimizing the interference power in the DD region makes decoding interference more 
difficult. This choice of action is not Pareto optimal. From the received power region representation in 
Fig. m we can see that the interference channel power should be maximized despite the values of the 
direct channel gain. It means that for each achievable desired channel power value, the interference 
channel power should be increased for easy interference decoding and removal. 

C. The Pareto boundary characterization 

From Thm. [5] and [31 we have presented the Pareto boundary characterization in ND and DD 
region. We can easily obtain the candidate set Q!^^ by reversing the role of Tx 1 and 2 from Vt^'^ in 
Thm. [21 Also, the candidate set f^""" is shown to be the following [16]: 



n^^ = {Wi, W2, Pl=P2 = Pma.} 

where Wi, W2 are defined in Thm. [21 



(24) 
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By definition in p9|) , the candidate set of the Pareto boundary B{TZ) is the union of the candidate 
sets in each decoding region. Hence, we have characterized the Pareto boundary B{TZ). In the 
candidate sets of the Pareto boundary, the beamforming vectors are parameterized with positive 
real scalars < Ai, A2 < 1- By varying Ai, A2 from zero to one and P2 from zero to Pmax, we obtain 
aU beamforming vectors that may attain the the boundary in each decoding region and in turn the 
overah Pareto boundary. Intuitively, it means that the boundary attaining beamforming vectors in 
each decoding region exist only in a two-dimensional subspace, spanned by the direct channel and 
the interference channel, in a A^-dimensional signal space. 

As a direct application of the Pareto boundary characterization, we characterize the maximum 
sum rate point in the following section. Since the maximum sum rate point is always on the Pareto 
boundary, the candidate set of the maximum sum rate point is therefore a subset of the candidate 
set derived above. We reduce the size of the candidate set by eliminating beamforming vectors in 
the candidate set that achieve a smaller sum rate than other vectors in the set. 

VII. The maximum sum rate point characterization 

In this section, we characterize the candidate sets of the maximum sum rate point by first 
illustrating that Txs. should always transmit with full power, in Section IVII-AI Then, we study 
the candidates sets of the boundaries of ND and DD regions by eliminating vectors that attain a 
smaller sum rate than other vectors in the candidate sets and obtain the candidate sets of maximum 
sum rate point in B{TZ"'^) and BiJZ'^'^) respectively, in Section [VII-BI and IVII-CJI 

A. Full power transmission 

We observe that the maximum sum rate point is attained by maximum transmit power at each 
transmitter. To see this, we combine the power constraints and beamformer norm constraints: 

||wif < Pi. (25) 

Assume that the sum rate optimal beamformer is not transmitting at maximum power: ||wi|p = 
p < Pi- We can choose a beamformer w^ = w^ + ee^'^Yljiha where e is chosen such that ||w^|p = Pi 
and (j) = arg(h^Wi). Notice that |h^w^p > |h^Wip and |hj|w^p = |hj^Wip. Or, we can choose 
w^' = Wi + e'e-''^ n^hji with (j)' = arg(h^Wi) to increase |h|[wip and keep |h^Wip constant. Thus, 
it contradicts that w^ is on the Pareto boundary. From now on, we set Pi = Pmaxji = 1,2. Note 
that the argument above is limited to non-parallel channels, for parallel channels (e.g. hj^h^i = 0), 
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it reduces to SISO-IC where the maximum sum rate point is attained by one Tx transmitting with 
full power whereas the other Txs. transmit at less than full power [13]. 

In the following sections, we characterize the candidate sets that attain the maximum sum rate 
point. Note that, the candidate sets attaining the maximum sum rate point is a strict subset of those 
attaining the Pareto boundary. The sum rate metric does not distinguish between Tx 1 and 2's rate 
and therefore we can identify a much smaller candidate set as illustrated in the following sections. 
This is particularly useful for system optimization which does not put emphasis on user fairness. 

Note that the computation of the global optimal solution 

ui* = argmax^"^(wi,W2) (26) 

Wl,W2€5 

is in general NP-hard [21], even though there exist channels for which the solution is easily obtained 
(e.g. orthogonal channels). Here, we would like to reduce the search space and characterize the 
solutions set. 

B. The maximum sum rate point characterization in the ND region 

Theorem J^: The candidate set of maximum sum rate BJ^'^ denoted as ^2"*^, hence u* C ^2'"^ C ri"'^, 
is given by 

^^^ = {Wi, W2, P^ax, Pma.\ (27) 

where r^"** is the candidate set of Pareto boundary B{TZ"'^) in (|2(J|) . In particular, VVi is the following 

set with cardinality three: 

Wi = |^,^,w,(Af))| (28) 

L||nii|| ||n2i|| J 

with Ai = —jrr — MO ^ '^■\ H^J I , — iri — ir>- The candidate set VVo is a set of beamforming vectors 

characterized by a parameter A2 in a smaller range than the range in W2: 



ni2h22 , r. r- n5L2h2 

in h II + V 1 - -^2,1^, 

|lll2n22|| ||lli2'l22| 



W2 = <! W2 e 5 : W2 = VA^^^i^ + yr^^Ji^ ; Af < A2 < Xr (29) 



where A^^^ = J ^mih I ^^ ^ parameter that gives the beamforming solution towards channel h22 
and W2(A2 ) = ,_ . v^ + ^ °. Vb for some eigenvectors Vq, v^ and positive scalars a, b. The vectors 
Vq, Vf, are the most and least dominant eigenvectors of the matrix S = h22h|^ — ^2ih^2hf^. 

Proof: See Appendix IXI-DI ■ 

Remark 4-' Note that for some channel realizations and chosen W2, wi(Ai ) may be equal to 

the maximum ratio transmission solutions m!^",, or j^^. But we distinguish between them in the 

||nii|| ||h2i|| " 

candidate sets to illustrate that for most channel realizations and W2, wi(A]^ ) ^ \ iTiriT' iIiriT f ■ 
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IhfaWspF; 



'Wi(Af^) 



^ "11 




f] 



nd 



{Pi = Pr 



'max,Wi G Wl,W2 G W2,0 < F2 < FmaJ 



|h22W2| P2 



Fig. 5: The illustration of the candidate set of the maximum sum rate point of ND region in dark 
grey and the candidate set of the Pareto boundary B{TZ^'^) in light grey. The cardinality of the 
candidate set for wi of the maximum sum rate point is only three, conditioned on W2. 



It is interesting to see that the sum rate optimal beamforming vector of Tx 1 is either the 
beamforming vector towards the desired channel or the beamforming vector towards the interference 
channel or a beamforming vector that balances the interference decoding rate Ti and the treating 
interference as noise rate Di in a weighted manner with weights ci , C2 which depend on the choice 
of the beamforming vector at Tx 2. 

Comparing the candidate set of the Pareto boundary to the candidate set of the maximum sum 
rate point of the ND region, namely Q"'^ and 17"'^, we observe the following: 

• For each Tx i, the candidate set Q"'^ consists of only three closed- form beamforming vectors 
whereas il"*^ consists of a set of beamforming vectors characterized by a real- valued parameter 
spanned between zero and one, as shown in Fig. \5\ 

• An interesting question rises: what are the conditions of each of these potential sum rate optimal 
solutions being sum rate optimal? We give the discussion in Section [Villi 
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C. The maximum sum, rate point characterization in the DD region 

In this section, we compute the candidate set that attains the maximum sum rate point of the 
DD region. 

Theorem 5: The candidate set of the maximum sum rate in R'^'^ is 



n" 



\^1 T^2 ' ^maxj J~^max f l'^"/ 



where for user i = 1,2, the sum rate optimal beamforming vectors are either a hnear combination 
of two orthogonal vectors or maximizing the desired channel power or a specific vector: 

^''=l^-^'^'(^^)| (31) 

^^ = 1^^ ■ ^WITT + v/^^^^^StTT' ^^ ^ ^^ ^ ^r4 (32) 

where Af^^ = n '^|!;^'' and \f = —-^, ,|'"t^''", » , / ■ 

Proof: see Appendix IXI-EI ■ 

Remark 5: Note that Wi(A^) may not be a element of Vi because A^ may not be smaller than 
A^*"* and in this case Vi is empty. The vector Wi(A^) is a beamforming vector that balance the 
interference decoding rate Ti and the treating interference as noise rate Di in a weighted manner. 
See Appendix IXI-EI for more details. 

In Fig. [Ul we illustrate the reduction of the candidate set of the maximum sum rate point of the 
DD region, in red, comparing to the candidate set of the Pareto boundary of the DD region, in blue. 
As shown in Fig. [6l the beamforming vectors in V achieve channel powers that are in the direction 
of minimizing the direct channel power while maximizing the interference channel power. 

To summarize, we obtain the candidate set of the maximum sum rate point in the ND and DD 
region, in Tlim. 3] and Thm. [5] respectively. We can exchange the role of Tx 1 and 2 in Thm. 3] to 
obtain the candidate set of maximum sum rate point in the DN region, Cl'^"'. For the NN region, 
the candidate set of the maximum sum rate point is identical to the candidate set of the Pareto 
boudary, il""". Thus, we can have candidate set of the maximum sum rate point of MISO-IC-IDC 
as Q: 

Sl = ^"^(Jfi'^"|JO'''*|Jl]™. (33) 

In the next section, we apply the results obtained above: from the candidate set of the maximum 
sum rate point in different decoding structures, we identify the conditions in which the MRT 
strategies are sum rate optimal. Such strategies are attractive because of their simplicity and the 
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|hf2W2|2F2 



W2(A^) 




llhii 



|hfiWi|2Pi |hgw2pP2 

^dd = {0 < Pi, P2 < P™a., Wi G Vl, W2 e V2} 
fi<i'^ = {Vi,V2} 

Fig. 6: The illustration of the candidate set of the maximum sum rate point of the DD region, in 
red, and the candidate set of the Pareto boundary B{TZ'^'^) in blue. If \f < Xf^'^, then the candidate 
set consists of the set Vi and w,; = jrcr^ as illstrated by the dark grey area in the figure. When 



||h,i 



A^ > Xf^^, the candidate set becomes 3 beamforming vectors: 



,w 



(Xf) and 



w,: 



MRT optimality conditions answer the following two interesting questions: When is selfishness sum 
rate optimal? When is interference amplification sum rate optimal? 



VIII. MRT OPTIMALITY CONDITIONS 

In this section, we investigate the conditions in which the MRT strategies at both Tx 1 and 2 are 
sum rate optimal. For clarification, MRT strategies can mean two strategies, one to beamform to 
the direct channel ha and the other to beamform to the interference channel hji. We characterize 
the MRT optimality conditions in terms of the separation between the desired channel and the 
interference channel, Of. 



cos 



I j^ 1 1 1 1 i^' 



(34) 



Theorem 6: The MRT optimality conditions for decoding structure R are: 



I 11 . " ii , 1 1. ^^ 11 I is optimal if and only if 

V ll"ll|l Il'l22|| ) '^ •' 



cilln^ihnf (l + ||h22fP™a.)||hrii2 



111 



< COS (^i) < ,, , 11^ no ,,^ ^^ Miu iic> (35) 



C2||h2iP-2VSIZi|hghn| + ci||hnP ' '" (1 + ||hi2Pcos2(02)P„.ax)||h2iP 

• ( llhT'' TlhT') ^^ optimal if and only if 

lln Vi I|2 

11, 112^/1 , ii, ii2o ^ J-i2inii , , 

h2i < 1+ h22 Pmax) ^ , n, n, TTTTT, 36 

1 + ||hi2rcOs2(6l2)Fmaa; 

^1^^^^ ^1 = ||h,.||^cor"(e;)p„.„.+i and C2 = ||h,,pHl.+i - 

Proof: See Appendix IXI-FI ■ 

Now, we provide the MRT optimality conditions for R*^*^. 

Theorem 7: The MRT optimality conditions for R*^** are: 

Wi = -!^ is optimal if -^ - 1 < ||hiif cos2(0,) < j\hijL (37) 

ll%i|| 5jj l + 5ji 

h 

Wi = " is optimal if (1 + (7jj)||hii|| < ||hji|| cos {6i) (38) 



for i,j = 1, 2 and gij = |h^Wjp. To be more specific 



(l+||h22f)||hi 




cos^(0i)>^'^" ''"'"^'' 



llVi 112 

is optimal if and only Hi „ „ (39) 

Cos2(02)>^" ''" ^" ''" 



V. 112 



is optimal if and only if 

h22f COS2(02) < llhnf ||h22f COS^ei)cOs2(02) < H ^21 IP 11^22 f COS^ (0; 



l + ||h22PcOs2(02) 

hi2n|hnfcos2(ei) 



(40) 



|h2ir-||hnrcos^(0i)< ||hnr||h22rcos^(ei)cos2(02) < , ,,^ ,,2 2/^ 

1 + ||nii|p cos^(6' 



Proof: see Appendix IXI-GI ■ 

IX. Simulation Results 

In this section, we provide simulation results for the proposed parameterization. By varying the 
beamforming vectors and power allocation, according to the proposed parameterization, we plot the 
achievable rate region for each decoding structure for a particular channel realization in Section 
IIX-AI The maximum sum rate point and the MRT points in each decoding structure are plotted on 
the corresponding achievable rate region. In Section IIX-BI we compute the empirical frequency of 
MRT strategies in R"*^ and R'^'^ averaged over 500 channel realizations. In Section [IX-Cl we allow 
the channels to be correlated and we see that the sum rate optimal decoding structure changes with 
the strength of the interference channel, agreeing with the observations for SISO-IC. 
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A. Achievable rate region and maximum sum rate point 

In Fig. [71 we plot the achievable rate region of the decoding structure R"'^, R'^", R'^'' and R"" in 

Fig. [Tal I7b| [7clandl7d | respectively. We assume A^ = 3 transmit antennas and SNR=OdB. We exhaust 

Ai and A2 to take 20 values between zero and one, inclusively. For each pair of (Ai, A2), beamforming 

vectors wi,W2 are generated and the corresponding rates with transmit powers Pi,P2 are plotted. 

Depending on the candidate set in each decoding structure, the transmit powers can be less than 

maximum power or equal to the maximum power Pmax- For example, in R"", maximum power is 

used: Pi = P2 = Pmax. whereas in R"^ Pi = P^„, and < P2 < P^ax and in R'''*, < Pi,P2 < 

Pmax- For simulation purposes, we allow the transmit powers to take 10 values between and Pmax, 

inclusively. The rate points plotted are achieved by the proposed Pareto boundary parameterization 

and the red asterisk is the maximum sum rate point by employing the maximum sum rate point 

parameterization where as the red square is the MRT strategies: ( wi = .n^^^, , W2 = n^^^ii ) in R"*^; 

fwi = |TTrTT,W2 = iTTri^l in R'*" and fwi = wrh,^2 = ttttt-I in ^'^'^ and R"". 
V ^ l|hii||' ^ I|hi2||y V l|hii|r ^ I|h22||y 

B. Empirical Frequency of MRT strategies 

In Fig. [HI we demonstrate the variation of the empirical frequency of MRT strategies in R"'^ when 
SNR increases. In R"**, the empirical frequency of beamforming vectors pair ( iThT-, \\\^^\\ ) i^ 50% 
when SNR goes to infinity where as the empirical frequency of (]Th^, TTirir) is 10%. It shows that 
with R"'', in high SNR, Tx 1 should amplify interference signal by beamforming at the interference 
channel and Tx 2 should amplify desired signal by beamforming at the desired channel and by doing 
so, it achieves maximum sum rate on average 50% of the channel realizations. 

In Fig. [HI we compare the maximum sum rate and the rates achieved by MRT strategies in R'^'^: 
('^i = PIT1I''^2 = pf]|) and (wi = p^,W2 = ptfl|)- In the x-axis, we plot the percentage of 
the average maximum sum rate whereas the y-axis is the percentage of channel realizations such 
that the MRT rates are less than a certain percentage of the maximum sum rate. Simulations show 
that the sum rate achieved by MRT strategies are less than 20% to 80% of the maximum sum rate. 
Since maximizing direct channel gain or interference gain do not reach maximum sum rate, it seems 
to imply that in R'^'^, interference should not be maximized or minimized and should be balanced 
instead. 

In Fig. [TUl we plotted the averaged difference between the maximum sum rate and MRT strategies 
in R"*^ and R'^'^ respectively. From Fig. llOal the rate difference between MRT strategies and the 
maximum sum rate decreases with SNR and reaches to about 2% and 4% at 40dB SNR. Thus, even 

20 



i.a 
1.6 
1.4 
1 2 


A 


RT 






'^ax. su 


m rat^ point 














\ 




















\ 










1.0 










































6 
















y 






4 
















\ 






0.2 
















\ 




















\ 






- 





2 


4 


6 


8 1.0 1.2 1.4 1.6 1 


8 





















^ 




^ 










































^^ 


^^ 






















^^ 












^ map(. surrt fate point 


































MRT 




















U 2 U 


4 b U 


8 1.0 1.2 1 


4 1 


fa 1. 



Rl 



Rl 



(a) Achievable rate region of R"'': proposed pa-(b) Achievable rate region of R''": proposed pa- 
rameterization achieves the Pareto Boundary andrameterization achieves the Pareto Boundary and 
maximum sum rate point. maximum sum rate point. 





















S 


1.8 

1.6 

1.4 

1.2 

1.0 
(M 

°^0.8 
0.6 
0.4 
0.2 








































max. 


sum r^te po 


X 








^■^ 


^ 






























N> 


\, 










^^ 
























MF 


T 




\ 








































^ 
























































, 


















- 




















\ 




















Tl 

^ 


1RT 












\ 






















































\ 






0.2 


0.4 


0.6 


0.8 1.0 1.2 1.4 1 

Rl 


6 1 







2 


4 





6 


8 1 
D1 


1 


2 


1 


4 


1 


6 1 





(c) Achievable rate region of R'*'*: proposed pa-(d) Achievable rate region of R"": proposed parame- 
rameterization achieves the Pareto Boundary andterization achieves the Pareto Boundary and maxi- 
maximum sum rate point. mum sum rate point. 

Fig. 7: Achievable rate region of different decoding structures. 



if the empirical frequency is about 50% and 10%, MRT strategies only lose about 2% and 4% of 
the maximum sum rate of R"'''. On the other hand, from Fig. I10b| it shows that the rate difference 
between MRT strategies and the maximum sum rate in R*^*^ increases with SNR and we conclude 
that MRT strategies are not sum rate optimal in R'^'^. 
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Fig. 8: MRT optimality in R"*^ when SNR increases: interference should be maximized half of the 
time when SNR goes to infinity. 
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Fig. 9: Maximum sum rate plots in diflFerent channel realizations. MRT strategies are not sum rate 
optimal in R'^'*. 
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(a) The averaged sum rate diflterence between maximum sum rate point 
and MRT strategies in R"''. 
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(b) The averaged sum rate difference between maximum sum rate point 
and MRT strategies in R'*''. 

Fig. 10: The averaged sum rate difference between maximum sum rate point and MRT strategies in 
R"'' and R'^'^. 
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C. Correlated Channels and Sum Rate optimal decoding struetures 

In this section, we assume a symmetric channel [22] in which the direct channels, ha, are i.i.d 
complex gaussian vector channels. The interference channel hji has a projection angle 9i with the 
direct channel h^^: 

Ih-^h'-l = llh-'llllh-'ll cos('6i'l (41) 

Moreover, we define the signal to interference ratio SIR as 

SIR=|^. (42) 

||"jj|| 

In Fig. [TTl we compare the sum rate achieved by R"", R'^'^ and TDMA. When the strengthen 
of interference channel increases and consequently SIR decreases, there is a transition from R"" 
to TDMA to R'^'^: treating interference as noise is sum rate optimal in low interference regime 
and then time sharing should be performed and then decoding interference is sum rate optimal in 
high interference regime. When the angle between the interference channel increases to 9 = O.lSvr, 
about 27 degrees, there is a direct transition between treating interference as noise and decoding 
interference. Thus, when the direct channel and the interference channel are more apart, time sharing 
is not sum rate optimal and outperformed by the other decoding structures. 

In Fig. [121 we compare the maximum sum rate achieved in different decoding structures with 
TDMA when the system SNR increases. When the interference channel is as strong as the direct 
channel SIR= 1, treating interference as noise is sum rate optimal in all SNR range. When the 
interference channel power increases SIR""*^ = 5, 10,20, both Rxs. decoding interference is sum rate 
optimal in low SNR whereas one Rx treating interference as noise and one Rx decoding interference 
is sum rate optimal in high SNR. 

D. Performanee of suhoptimal algorithm 

For illustration purposes, we propose a very simple transmission strategy with only finite low 
number of beamforming vector choices. This transmission strategy is inspired by the parameterizaiton 
of each decoding structure. We propose to select only two beamforming vectors in each candidate 
set. Based on channel states information, we compare the sum rate performance of these eight 
beamforming vectors and choose the beamforming vector and the corresponding decoding structure 
which achieves the highest sum rate. 

. NN region: f JI^^, ^i^) and f^, ^a^V 

^ Vl|n2lhll|| ' ||ni2h22||/ Vll'^llll ll'^22||/ 
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Fig. 11: The averaged sum rate of different decoding structure comparing to TDMA in symmetric 
channel when the strengtli of the interference channel increases. 
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Fig. 12: The averaged sum rate of different decoding structure comparing to TDMA in symmetric 
channel when the system SNR increases. 
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ND region: (t^, t^) and (j^, j^ 

o Vl|h2l|| ' ||h22||/ V l|hii|| ' ||h2 

DN region: f^, ^) and f ^, ^ 

o Vll'^iill I|hi2||y V l|hii|| ' ||h2 



. DD region: f J^^ hi^) ( hi^ ha^) ( h2^ ha^) and f ^, ^1 • 
° Vllh2i|r ||hi2||y ' vilhiilr 11^2211/ ' Vlli^2i|r 11^2211/ vili^iilr lli'iall/ 

• TDM A: a time sharing scheme between single user points and w^ = ttjj^- 

In the NN region, we propose to choose either the interference nuhing solution or the desired 
channel gain maximizing solution. It has been shown in previous literature that in MISO-IC-SUD 
and low SNR regime, maximizing desired channel gain is sum rate optimal whereas in high SNR 
regime, interference nulling is sum rate optimal. 

It was shown in SISO-IC [3], [lOjthat the DD scheme is sum rate optimal among all four decoding 
structures when the strength of both interference channels are strong and the DN or ND scheme 
is sum rate optimal when one interference channel is strong and the other interference channel is 
weak compared to the desired channel. The interference maximizing beamforming solution and the 
desired channel channel power beamforming solution are chosen in DN, ND and DD regions in the 
proposed algorithm to verify the analogy from SISO-IC to MISO-IC. 

In Fig. [131 we plotted the maximum sum rate achieved by different decoding structure and compare 
it with the proposed simple algorithm when SIR decreases. We see that when the interference is weak, 
it is sum rate optimal to treat interference as noise and when the interference strength increases, sum 
rate can be increased by allowing one of the Rx to decode interference and in the strong interference 
regime, both Rxs. decoding interference achieves the highest sum rate. Depending on the channel 
coefficients, TDMA may outperform R"" and R'*'^ in the medium interference regime. Note that the 
computation of the maximum sum rate point is NP-hard. However, we see the the proposed simple 
algorithm achieves nice sum rate performance with only five choices of beamforming vectors. 

X. Conclusion and future work 

The interference decoding capability brings additional freedom to the Rxs. which either decode in- 
terference or treat interference as noise. However, it is not trivial when Txs. should avoid interference 
and when should amplify interference power. To answer this question, we formulate the achievable 
rate region for a two-user MISO-IC-SUD. We provide an in-depth analysis of the achievable rate 
region, as a union of different decoding structures. The Pareto boundary is then characterized in 
terms of both power allocation and beamforming vectors. As a direct application of the Pareto 
boundary characterization, we characterize the maximum sum rate points. The candidate set to the 
maximum sum rate point is a strict subset of the the candidate set of the Pareto boundary. With 
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Fig. 13: Sum rate optimal decoding structures when the strength of interference channel increases. 



the maximum sum rate characterization, we derive the MRT optimality conditions which describe 
the conditions in which simple MRT strategies are sum rate optimal. We conclude the paper by 
providing simulation results which shed some insights into the question " When is selfishness sum 
rate optimal?". Results show that MRT strategies have only 2% to 4% rate loss comparing to the 
maximum sum rate point in high SNR if one Rx decodes interference and the other treats interference 
as noise. On the other hand. MRT is not sum rate optimal if both Rxs. decode interference. In 
symmetric channels, there is a transition in decoding structure from treating interference as noise 
to TDMA to decoding interference at both Rxs. when the strength of interference increases. 

The achievable rate problem for the A'-user MISO-IC-IDC is not simple as the number of possible 
interference decoding order increases exponentially with the number of users. It is not easy to see 
which decoding order is better than the others and whether this decoding order in this decoding 
structure is better than other decoding structures. This extension to the K-user case is currently 
under preparation. 
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XL Appendix 
A. Proof of Thm.Ul 

In this sequel, we are going to prove that the Pareto boundary in NN region, DN region and DD 
region are attained by rank 1 matrices. 

Lemma 1: In the NN region, where Rx 1 and Rx 2 treat interference as noise, the Pareto boundary 
attaining transmit covariance matrices are rank one. 

Proof: See reference [4], [20]. ■ 

To prove that the Pareto optimal transmit covariance matrices in DN and DD regions are rank 
one, we follow and modify slightly the proof from [20]. To facilitate the discussion, we define the 
following received channel power region for Tx i, assuming transmit power being one, 

$f = {(hf S.h,„ hf.Sihj,) : Si G cS} . (43) 

For each received channel power gji = h^S^hji, i,j = 1, 2, there is a set /Cjj such that each utility 
Uk, k £ /C'j, is monotonically increasing with received channel power gji, 

Uk {gu , Qji , gtj , gjj ) <Ukigii, g'ji , gtj , gjj ) (44) 

if <7ji < g'ji- Similarly, there is a set /CJ, such that each utility Uk, k E /Cjj, is monotonically decreasing 
with received channel power gji. 

Lemm.a 2: For an arbitrary fixed Tx i, if there exist a received channel power gji, j = 1,2, such 
that the number of utilities that are monotonically increasing and decreasing with gji are both larger 
than zero and the number of such received channel power gji is not larger than one, e.g. 

A/; = {/.:|/Cy >0,|/Ct,| >0},[[A/;[[ <1, 1 = 1,2, (45) 

then the Pareto optimal transmit covariance matrices with respect to these utilities are rank one 
and attain the boundary of the corresponding received channel power regions $1 . 

Proof: We proceed by separating the cases where the some utilities are increasing or decreasing 
with the received channel power gji, for i,j = 1, 2. 



. If 






7^ and /Cjj = 0, then the Pareto optimal transmit covariance matrix S* attains the 
boundary of the received channel power region $®. It is because for each utility Ukign, gji,gij,gjj) 
such that k E /C!j, if gji is not on the boundary of $f, we can choose a transmit cova- 
riance matrix S* such that g'ji = h|[S*hji > gji which increases the value of the utility 
Uk{gii,gji,gij,gjj) < Ukigu, g'ji, gij, gjj)- Since /CJ^ = 0, no utilities are decreased by modifying 
the transmit covariance matrix from S to S*. 
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If 






7^ and /C'j =0, then the Pareto optimal transmit covariance matrix S* attains the 



boundary of the received channel power region $?. It is because for each utility Ukigu, gji,gij,gjj) 
such that k G /CJj, if gji is not on the boundary of $*, we can choose a transmit cova- 
riance matrix S* such that g^^ = h^S*hji < gji which increases the value of the utility 
Uk{9ii,9ji,9ij,gjj) < Uk{gii,gji,gij,gjj). Since /Cjj = 0, no utilities are decreased by modifying 
the transmit covariance matrix from S to S*. 

According to the assumption ()^5|) . for each Tx i, there exist at most one j such that O, ^ 
and 






7^ 0. In this case, the received channel power with transmit power, gjiPi has value 
between and ||hji|pPmQa; which can be achieved by setting gji to be on the boundary of $? 
and letting Pi vary from to Pmax- Hence, we can put the Pareto optimal transmit covariance 
matrix S* on the boundary of the received channel power region $| by allowing power control 

n < p < p 

^ _ -^ 2 _ -^ max- 

From [20, Lemma 3], the received channel powers on the boundary of the received channel power 
region ^f are attained by rank one transmit covariance matrices, which completes the proof. ■ 

In the following, we are going to show that the Pareto optimality problem in DN region and DD 
region respectively satisfy the assumption in Lemma [2] and therefore the Pareto optimal transmit 
covariance matrices are rank one. 

1) In the DN region: Rx 1 decodes interference and Rx 2 treats interference as noise and the 
Pareto boundary attaining transmit covariance matrices are rank one. To see this, we recall the 
achievable rates in R'^": 



(46) 



■"1(511,521,512,522) 


= log2(l + giiPi) 






^^2(511,521,512,522) 


/ / 
= min log2 1 + 
V V 


512-P2 \ , 

n P ^ 1 j ' ^°^2 

giiPi + i/ 


f^ j 522^^2 
V l+52li^l 


and therefore we have 








4i = {i}; 


^ii = {2}; 


4i = W; 


4i = {2}; 


ICI2 = {2}; 


/ci2 = W; 


^22 = {2}; 


)ci, = m. 



(47) 

Since the assumption (HSl) is satisfied, by Lemma El we have the Pareto optimal transmit covariance 
matrices in DN region as rank one. 
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2) In DD region: Rx 1 and 2 decode interference and the Pareto boundary attaining transmit 
covariance matrices are rank one. To see this, we recall the achivable rates in R'^'^: 

521 



ui{gii,g2i,gi2,g22) = min (logj (l + 511) ,log2 ( l + 
U2igii,g2i,gi2,g22) =min (log2(i + 522) ,log2 ( 1 + 



1+522 
ffl2 

1 + ffll 



(48) 



Hence we have 














4i = {i}; 


4i = {2}; 


4i 


= {1}; 


ici = {0} 




/Cl2 = {2}; 


^2 = {0}; 


^^22 


= {2}; 


42 = {1} 



(49) 

Since the assumption (|45|) is satisfied, by Lemma [21 we have the Pareto optimal transmit covariance 
matrices in DD region are rank one. 

B. Proof of Thm. 

The Pareto optimal beamforming vectors can be presented as the solutions of the following 
optimization problems, for some feasible value r^, i,j = 1,2, j ^ i, 

max Rj{-wi,W2,Pi,P2) 

Wl,W2,Pl,P2 

subject to Ri{y^i,^2,Pi,P2)>ri, ^ ^ 

(50) 

||wi|| = 1, ||W2|| = 1, 

< -Pi < Pmax , < -P2 < Pmax ■ 

With different decoding structure, the achievable rates Ri, R2 are substituted with different rate 
expressions. However, both the ND and DD case resemble closely to the maximization problem of 
channel power and the optimal solution can be characterized in a linear combination of some specific 
vectors. 

Lemma 3: Define a maximization problem of channel power in the following: 

max t 

subject to |u-^wp > ut, |v^wp > vt (51) 

||wf < 1 
for some arbitrary fixed scalars u > and v >0. The solutions w* must be in the following set W, 

W = |w G C^^i : w = Vr^T^ + yr^^jl^, <fi<l]. (52) 

I lliiuvii l|n^v|| J 

Corollary 1: Reversing the signs of the inequality in (|5ip to |u-^wp < ut, |v^wp < vt does not 
change the characterization of the solutions in (|52|) . 
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Proof: We proceed by writing the Lagrangians of the problem, with Lagrange muhiphers A 
[Ai, A2, A3], 



Lfw, X) =t 



A.( 



ut 



lu^wP 



A2( 



vt 



Iv^wP 



A. 



w 



11 



(53) 



Now, we compute the vanishing point of the Lagragian derivative which is a necessary condition of the 
optimal solution, — 1^^ ' = Aiuu'^w+A2vv^w— A3W = 0. We can write Aiuu^w+A2vv^w = A3W 
and adjusting the constant scaling, we have Ai||u|pT]^Ww + A2||v|P|St2W = A3W. Therefore, the 
eigenvector w is a composition of its projection on u and v: 



Ai||u|pn„w + A2||v|pn„w = A3W. 

Since A^ > 0, i = 1, 2, 3, we can write, for some complex-valued /ii, /i2, 

r^^ Hull r*^ ||v|| 

w = 



/ilFT+M2rT 



Ml 



Hull 



+/i2 7 



Now we define the set of beamforming vectors that satisfy (|55|) ,U = < w : w 

Then, we show in the following that U is a subset of W in Lemma [3 We start with 

u V 

U V 



(54) 
(55) 

,/il,/i2 G C 



(a) ^1 



|U|| 

' ^J'l 



n„ 



1U.112 



M2- 



V^U 



■e-J''^'"'n„u 



where (a) is due to v 



Hvullvll 



lTTTT 



|u|| 

nt,u 

||nt,u| 

and (puv 



+ ^21 



re 



n„u + n,iu 



(56) 



Zl 



Z2-: 



\n^u\\ 



arg(v^u). The parameter zi is complex and Z2 is real; 

(a) 

the values of zi and Z2 are scaled such that ||w|| = 1. Notice that |u"^w| = |zi||n,i;u|| + 2;2||ni^u||| < 



|zi|||n„u|| +Z2||n^u|| and |v"^w| 



zi- 



|n„u|| 



\zi\ 



|n„u|| 



Note that equality at (a) when zi is real 



and the phase of z^ does not affect |v^w|. Hence, Zi can be chosen real and since the two basis are 
orthogonal, the power constraint of w is satisfied when z\ + c^ = 1 and thus, we can write zi = yjji 



and C2 = \J1 — /i for < ;U < 1. 

If the signs of inequalities are reversed and we have |u-^wp < ut, |v^wp < vt in the constraints 
in (fSTj) , then the corresponding signs changes to minus from positive in the Lagrangian but does not 
affect the discussion above and the characterization of the Pareto optimal solutions holds. ■ 



31 



Substitute the rate definitions ([3]) into (j50|) and maximizing Ri subjecting to a constraint on R2, 
we let Z22 = |h^W2pP2 and z 



12 ~ |hf^W2pP2 and we focus on the subproblem that concerns wi, 
t 



max 

"Wl 



subject to 



\hgwi\^>{zl2 + l)t, 

|hfiWi|2>(z22 + l)t, 



(57) 



l|wi|| < 1- 
This has the same formulation as in Lemma [3j By substituting u = h2i and v = hn, we obtain the 

characterization of the Pareto optimal beamforming vectors as a linear combination of the vectors 

n2ihii and H^j^hn. Now we reverse the optimization order : maximize R2 subject to a constraint 

on Ri. After some manipulations, we obtain for some zfi, Z21, 

|hiw2pP2 



max 

W2,P2 



subject to |h22W2pP2 ^ 

|hf2W2pP2 < 



^21 



2'-i - 1 



1, 



(58) 



^11 



2'-i - 1 

IW2II < l,0<P2<i^rr 



By Corollary [11 we have 



for < ;U2 < 1. 



W2 



Hioh 



12"22 



Bioh 



+ VT 



nfnh 



M21 



12"22 



12"22 



nfoh 



(59) 



12 "22 1 



C. Proof of Theorem\^ 

The following proof is similar to the approach in Appendix IXI-BI to avoid repetitions we only 
highlight the main differences in the following. The Pareto optimality problem in the DD region is 
written as a maximization of Ri subject to R2 > ^2 and after some manipulation, we focus on the 
subproblem optimizing wi: 

max t 

Wl,Pl 



subject to |h^wipPi > t, 

IhfiWipPi < 



212-P2 



-1, 



2r2 _ 1 

\K^l\^Pl>t{zl2 + l), 

II Will <l,0<Pi<P^a. 



(60) 
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for some z^2 j -^22 • Similar to Lemma [3l we write the Lagrangian and set the derivative with respect 
to wf^ to zero, and obtian {Pi{Xi — A2)hiih|^ + A3Pih2ih^) wi = A4W1. We add PiA2||hiipwi to 
both sides and obtain 

(PiAillhnfnn + PiA2||hnf Hj^i + AgPiUhsif n2i) wi = (A4 + PiA2||hnf )wi. (61) 

Hence, we see that the optimal solution is composed of its projection on the subspace spanned by 
hii,h2i and the orthogonal subspace of hn, which can be represented by the following: 

Hiihsi n5Lh2i 

"^^^^ Min h II +^^ iimh II ^^^) 

||-lliin2i|| P-liin2i|| 

for some /ii, /i2 S C and |/iip + |^2p = 1- 

Similar to the arguments before in Appendix IXI-Bl we omit the details here to avoid repetitions. 
The values /ii, /i2 can be chosen real- valued and the Pareto boundary of R'^'^ attaining beamforming 
vectors are 

wi G Vi = I w = V^ "''!""' II + ^^^^ iini'l!"' II , < ^1 < 1 1 ■ (63) 



D. Proof of Theorem [^ 

Before we go into details of computing the candidate set of the maximum sum rate in the ND 
region, we present the following lemma which holds importance in the discussions later. 

Lemma 4-' Consider two functions fi{x), f2{x) where fi{x) is concave and f2{x) is linearly incre- 
asing with X G Af C M. Define 

x*i = argmax/i(x), X2 = argmax/2(a;) and X = {x G ^ : /i(x) = /2(a;)}. (64) 

If X = argmaxmin(/i(a;), /2(x)), then 

xC{xl,x*2,X} (65) 

Proof: Notice that since fi{x) is concave, there are at most two intersection points. If there 
are no intersection points, then x C {a:i,.T2}. If there is one intersection point x, then one of the 
following orderings is true: 

• X < xl < X2 '■ since there is only one intersection point, /2(a^2) < /i(^2) ^^'^ thus x = x\. 

• x\ < X < x\ : \ei x~ < X < x'^ and we have fi{x~) > fi{x) > /i(x^) and f2{x~) < f2{x) < 
f2{x~^)- Thus, X = X. 

Note that f2{x) is linearly increasing with x and therefore X2 is the boundary of X and thus X2 > x^ 
and X2 > X. If there are two intersection points, then we have x C X. ■ 
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Note that the global optimal solution u* must be in the solutions set of B(R"'^): 

u* C n""^ (66) 

where Q,^'^ is defined in Thm. [21 Thus, we can refine the constraint set in the maximum sum rate 
problem in the ND region to 

^"'*(wi,W2)= max C2(w2) +min{ri(wi,W2),Z?i(wi,W2)} (67) 

(wi,w2)en'"' 

which can be decomposed to the following: 

max <C2(w2)+ max niin{ri(wi, W2), -Di(wi, W2)} > (68) 

w2eW2 [_ wieWi J 

where the inner maximin problem is maximized over wi for each v/2 ■ With each given W2, we 

define Ci = tt-h — ffrf rr and C2 = ttr — i^rr rr and we rewrite Ti and D-i to the following: 

Di = 2^^ — 1 = ci|h{^wip and Ti = 2^^ — 1 = C2|h2iWip. Thus, the maximin problem in ([68l) is 

equivalent to 

max min{ci|h^^wip, C2|h2]^wip}, for arbitrary fixed W2 G S. (69) 

W16W1 

Lemma 5: The candidate set of the maximin problem in (|69l) and therefore maximum sum rate 
problem in (f68|) can be reduced from Wi G il"'^ to VVi which is a candidate set with cardinality 
three containing at least the maximum sum rate solution: |VVi| = 3, 

Wi = |^,^,w,(Af))| (70) 

L||hii|| ||h2i|| J 

with 

Xib) ^ ^IP^lhlllP .7^^ 

' C2||h2i||2-2v/^IZJ|hfihn|+ci||hn|P' ^ ' 

Proof: Because of the formulation of wi € Wi, we can write the beamforming vector as a 

function of a real-valued parameter Ai, wi(Ai). Using the result in Lemma [3l we define /i(Ai) = 

ci|hf]^wi(Ai)p and /2(Ai) = C2|h2iWi(Ai)p. It is easy to see that /i(Ai) is concave in Ai and 

/2(Ai) is linearly increasing with Ai. The function /i(Ai) = Ci|hf(wi(Ai)p attains maximum when 

''^i(''^i) = |n^"|| ■ Similarly, f2{^2) attains maximum when u'i(Ai) = ttj^^- Now we compute A-^ 

which satisfies ci|hfiWi(Af^)|2 = c2|hgwi(Af^)|2. 



To proceed, we compute the channel powers |h^wip = I y '^i ||n2ihii|| + y \ — X^ '||n^-^h 



Ab) 



kW|lTT..U.JI^./l _ \W 



111 



and |h2iWip = A-^ ||h2ip. Notice that ||n2ihii|| = ||h2i|| cos((/)) and ||n2ihii|| = ||h2i|| sin((/)) where 
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Ih^ihiil = ||h2i||||hii|| cos((/)) and by definition cos((/)) is positive. Rewrite A-^ = cos^{6) where 
< 9 < 7r/2. Thus, we can rewrite the channel powers to 

|hfiwi|2 = ||hiifcos2(0 -(/)), 

(72) 

|hfiWi|2 = COs2(e)||h2lf. 

Thus, ci|hfiWi(Ai )p = C2|h2iWi(Ai )p is equivalent to y^\\hii\\cos{9 — (/)) = ,,/c^cos(6')||h2 



^21 II 

which is due to the fact that cos(0 — (/)) and cos(0) are by definition positive. Putting the sinusoids 

one side and we obtain '^°pg|-gf = ^^ L ^ | . If we expand cos{6 — (/)) = cos{9) cos{(f)) + sin(0) sin((/)), 

we have 

tan(0) = V^l|h2i||-Vcr||hn||cos(0)_ ^^^^ 

v^||hii||sm((/)) 

Use the Pythagorus theorem, if tan{9) = | then cos{9) = . ^ and therefore 

, . ^. IITT-L U II 2 

'' = ■=" "" = e.||h..|P-2,/ere;|hgh„|+e.||h„|P ^ <'"' 

■ 

If we reverse the maximization order in ([67|) . we obtain: 

i?" (wi,W2) = max max min{C2(w2) + Ti(wi,W2), C2(w2) + I?i(wi, W2)} . (75) 

wievvi w2eW2 

Lemma 6: The optimal solutions to R^'^ in (|75p can be reduced from W2 G f^"'^ to VV2, a set of 
beamforming vectors that includes the beamforming vector towards the desired channel h22, 



IIl2h22 r ^ 11^2^2 

in h II + V 1 - ^^-iifTTTT , 

|lll2n22|| ||lli2'l22| 



W2 = W2 G 5 : W2 = v^^^i^ + v/T^A^^^i^ ; A(^) < A2 < A-- (76) 



where A^^'^ = ,,} ^^mP , is parameter that gives the beamforming solution towards channel h22 and 

^ ||hi2||||h22|| 100 ^^ 

W2(A^0= / - Va+ ^- ^ Vb (77) 

Va + ^ V ffl + fo 

for some eigenvectors Va,Vf, and positive scalars a,b. The vectors Va,Vf, are the most and least 

dominant eigenvectors of the matrix S = h22h|^ — ^^h.i2^r2- 
Proof: R^'^ in (|75|) is equivalent to compute 

max min (1 + 521 + IhiwspPs, (1 + \h^22^ 2? P^a.) U + , , |,/''|2d ) I ' (^8) 

W2|wi L V l + \n.';2^2VPmaxJ ) 

where the notation maXvv2|wi denotes maximization over W2 for some given wi and can be decom- 
posed into the following two subproblems: 

max^,|„, 1 + 521 + \^22^2?Praax if \^22^2?Praax > §77 " 1 + ^f^\h.^2^2? Praax 

max^,l^,(l + |hf2W2pP2) (1 + i^ih^r^^pp ) if IhgwspP^a. < ft " 1 + fl\^?2^^? P^ax ' 

(79) 
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The first subproblem has optimum solution i n ^^m . For the second subproblem, the optimal A2 must 
be in the region A = {A2 : A^^^ < A2 < A^^} where W2(A^'^^) = ^ and |hgw2(Af^)pP^^^ = 
^^ — 1 + ^^|hf^W2(A2 )pPmaa;- To See this, we write the metric in the second subproblem as a 
function of A2: 

511 



F(A2 



l + |h^2W2(A2)rP™ax) 1 + 



l + |hf2W2(A2)PP„ 

Assume A"*", A~ ^ A, in particular, A"*" > A2^^^ and A~ < Ag , we have 

|hf2W2(A+)|2 = A+||hi2f >A"^^||hi2|P 

|hiw2(A+)p < |hiw2(A-^^)p = llh. 



i22| 



BO) 

Bl) 



(6) 



Thus, any A^ achieves -F(A~'~) smaller than F(A^^^^). Note that for any A~ < Aj , W2(A^) is not in 
the constraint set of the second optimization problem. It is because |h^W2(A)p is concave in A and 
attains the maximum at A2^^^ and |h|^W2(A)p is linearly increasing with A. Since A2 < Aj™^, for 
points A" < X^i^ we have |hf2W2(A-)|2 > 221 _ l + 22i|hf2W2(A-)|2. 

Unfortunately, the direct computation of A2 is tedious and does not give much insight. Here, 
we provide a cleaner method of computing W2(A('')) directly. Denote g = ^^. We compute the 
beamforming vector W2 such that 

|hiw2pF2 = 5 - 1 + 5|hf2W2|2p2. (82) 

Define S = h22hg - ffhi2hi2, S = S - {g - 1)1. 

wf Sw2 = (83) 

is a necessary condition for satisfying (|82|) . 

From the definition of S, we know that S is rank two with one positive eigenvalue and one negative 
eigenvalue [20]. Denote the non-zero eigenvalues of S by a and —b where a, 6 > 0. Employ eigenvalue 
decomposition on S and we have 

a 0ix(Af-2) 

S = [VaVb|V] -b Oix(Ar_i) [VaVb|V 







(N-2)x{N-2) 



H 



4) 
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where V is the A^ by A^ — 2 matrix with column vectors of eigenvectors of S that are orthogonal to 
Va, Vfe. With the same eigenvectors, we can write S as the following: 

a -(5-1) Oix(7V-2) 

S = [V,V6|V] -b-{g-l) Oix(W-l) [VaVblV 

-(5 - l)I(Ar-2)x(Af-2) 

Let a = a — {g — l),b = b — (g ~ I). The beamforming vector W2 of the following form 



H 



5) 



1 



Vi^2 



;V6 



6) 



Vffl Vb 

satisfies w^Sw2 = 0, where j = \/— T and (f) is a phase angle between to vr. Notice that if W2 has 

any power on the remaining orthogonal subspace spanned by V, W2 also satisfies (|82p but the value of 

|h|^W2p is smaller and therefore cannot achieve the maximum sum rate. It is easy to see the result by 

a 



direct computation: w|'^Sw2 = I -j^^a + 



e-Jf^M 



VaVfe 







^^-^'']'' 1 71^- 



pj4> _ 



rVfo 







for given angle (j). The formulation in (j86j) gives a family of beamforming vectors, each with a 
different value of (f). To fine the unique (j) and therefore W2 that maximizes sum rate, we rewrite the 
optimization problem in (|69p to 



1 , e-'* 

Va ^1, 



max |hgw2((/))p 

such that W2 ■ 
Define the following phase angles, (pa = arg(h|^Va), <?l)fo = arg(h|^V5) and therefore 



(87) 



//mfh^Vji 



+ 



m = a,b. 



otherwise. 

The optimization problem in (f87|l is therefore equivalent to 

1 



max |h^W2p = max 



1 e^'t' 

s/a Vb 



J 0b 



max 



Ihiva 



oj{<t>+4>b-4>a) 



Vb 



|hf2V6| 



Thus, the optimal phase angle cp is 



(88) 



To satisfy the norm constraint, we can scale the beamforming vector with a positive scalar, which 
does not change the direction of the vector. 

b 



W2 



;Va + 



eJ'I'd. 



V a + b ya + b 



^^b- 



9) 
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Note that if A2 > A™*"*, then yV'2 has only one element, i.e. -inj^. Lemma [5] gives the optimal 
candidate set for wi for arbitrary fixed W2 whereas Lemma O gives the optimal candidate set for W2 
for arbitrary fixed wi. Combining both Lemmas, we obtain the maximum sum rate candidate sets 
for both wi and W2. 

Remark 6: The authors in [23] provided a general solution of (|69|) . in the context of a multicast 
SNR balancing problem. The authors transformed the channel powers balancing problem to a 
weighted sum channel powers maximization problem for some positive weights Wi^W2- The optimal 
beamforming vector is then characterized as a dominant eigenvector of some matrices, depending on 
Wi,W2- However, the computation of such weights wi,W2 is not provided or trivial. In this paper, 
due to the beamforming vectors parameterization proposed in Thm.[2]and[31 we obtained the closed 
form solution of such channel powers balancing beamforming vectors. 

E. Proof of Thm. [3| 

In this section, we provide the proof of Thm. El We start by identifying four constraint sets of 
beamforming vectors 57*^°, 17^-*^, H^'', H^^ where the sum rate function is a different function in each 
set. In other words, when the beamforming vectors vary, the sum rate function being a sum of two 
minimum of rate functions, may change from one rate expression to another rate expression. The 
constraint sets are the set of beamforming vectors for which the sum rate function remains at one 
rate function. In the following, we provide the analysis for W2 but the sum rate function R'^'^ is 
symmetric with Tx 1 and 2. Therefore, we can exchange the role of Tx 1 and 2 and obtain the 
candidate sets for wi. The sum rate in the DD region is: 

R'^'^ = niinjCi, Ti} + min{C2 + T2} 

= min{Ci + C2, Ci + T2, Ti + C2, Ti + T2}. 

Z\ Z2 Z:^ Z4, 

We analyze each term and define Zi = 2^' , i = 1 , . . . , 4 

Zi = (1 + |hfiWi|2p)(l + |hf2W2|2p) 
Z2 = l + |hfiWi|2F+|hf2W2pF 

Z3 = 1 + IhfiWipF + |h|^W2pF 

^ _ / M,^,?P \ / , |hf2W2pP \ 

I 1 + |hiw2pP j y^ 1 + IhfiWiPP j 
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(90) 



(91) 



Notice that Z\ < Z2 is equivalent to Z3 < Z4 and Zi < Z3 is equivalent to Z2 < Z4. We summarize 
to the following lemma. 

Lemma 7: Let gu = |hf^wipPi and (/21 = |h2iWipPi. 

Zi ^ Z2 '^ Z^ < Z4 

^ (l + 5ll)|hf2W2|2<|hf2W2|2 
Zl ^ Z^ 44> ^2 ^ ^4 

^ |hiw2pP2<— -1 

511 
Proof: It is by direct manipulation of the definitions. ■ 

To facilitate representation, we denote the following two indicators and the corresponding candidate 

sets: 

^ ^ ' 1 if(l + gii)|hiw2pF2<|hf2W2pF2 
otherwise. 

, 1 if |hgw2p-P2 < ^ - 1 
otherwise, 

f^'^^ = {w2 : IIW2II = l,A = a,5 = fe} (94) 

Notice that Q!^^ , a,b = 0, 1, gives four constraint sets. By Lem.[71 we can decompose the optimization 
problem in (|90p to the following: 

R = log2 I max max nim{Zi, Z2, Z^, Z4} 

log2 fmaxwi max^^en" ZA ifA=l,i? = l 
log2 fmaxwi maxwjgnoi ^2) if A = 0, i? = 1 
log2 fmaxwi maXwaenio Z3] if A = 1, i? = 
log2 fmaxwi max^vjenoo ZA if A = 0, B = 
Now, we proceed with the proof of Thm. [S]in two parts: first, in Section [XI-E. 11 we identify the 
candidate sets for each of the subproblems Zi,i = 1, ... ,4; second, in Section [XI-E.21 we combine 
these candidate sets to one superset by eliminating beamforming vectors, which are on the boundary 
of the constraint sets, if they achieve smaller sum rate than other beamforming vectors. 

1) The candidate sets of subproblem,s Zi: The candidate sets for each of the subproblem is as 
follows: 

, If w = arg max^^gnii ^1 , then w G J7^ IJ f^^ U pfiT 
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. If w = argmax^^gnoi ^2 , then w £ W2(A^) U w^^ [j j^^. 

, If w = argmax^,gf2io ^3, then w G W2(A^) IJ ^^^ . 

• If w = argniaXw2gnoo Z4, then w G il"^ U VV|'^ 
where the constraints ^7 and Q^ are the set of beamforming vectors that satisfy the constraints 
by equality n^ = {w2 : (1 + gn)\hgw2\^ = IhfaWap} and f7^ = {w2 : |hgw2pP2 = |f^ - l}. The 
beamforming vector W2(A^) is the beamforming vector in Q^ that maximizes the desired channel 
power, W2(A2) = argmaXvv2gf2-4(l + 5ii)|h^W2p. The beamforming vector Wj is a unique vector 
that is a member of both ^7 and f7^ , W2 = ^ f] ^^ ■ Lastly, we have V2 which is a subset of V2 , 

V„ — /lir^ • , A~ n22hl2 1 n \~ n33hl2 \A ^ \ ^ \MRt\ ,„>,oro \MRT _ \il-22^'i-2\' 

V2 - |W2 . VA2||n22hi2|| + V i - A2pjp^, A2 < A2 < A2 I where A2 - pi;^^pr7^- 

Notice that beamforming vectors in 17^,^2^ satisfy the constraints in ([92|) with equality. To see 
this, we have the following observations: 

1. Zx is monotonically increasing with |h^W2p. If the constraints are not active, the optimal 
solution is m. ^^m . If the constraints are active, Zi is maximized over constraint set Q}^ and therefore 

||h22|| ' ^ 

the optimal solutions are in Q.^^^OP . 

2. 7,2 is monotonically increasing with |h{2W2pin constraint set ^^ . If the constraints are not 
active, the optimal solution is |Tj7f]T- There are an upper bound on |hf2W2p and an upper bound on 
|h|^W2p which in turn upper bound |h^W2p. If the constraint for |h^W2p is active, the solution is 
W2(A2 )• If both are active, the solution is Wg ■ Thus, the candidate set is ■^^V2(A2 ), W2 , itj^^ \ ■ 

3. Z3 is monotonically increasing with |h22'W2p in constraint set 17^*^. There is only one constraint 
that upper bound the value of |h^W2p, the optimal solution is in set 57 which at the same time 
maximizes |h^W2p, denote as W2(A2 )■ If the constraints are not active, we have ttjtt-- 

4. Z4 is monotonically increasing with |hf^W2p and decreasing with |h^W2p in constraint set 
V^^ . If the constraints are active, the optimal solutions are in Vt^ . If the constraints are not active, 
the optimal solutions are in V2- Similar to the case of ^"'^, any A ^ (A2 < A2 < A2™^) cannot 
attain maximum value of Z4 in f7°°. Notice that for any A > A'"''*, |hf^W2(A)p < |hf^W2(A^^^)p 
and |hgw2(A)|2 = A||hi2|p > A^^^||hi2|p = |hf^W2(A^''-*)|2 . Thus, Za^{\) < ZiiXf""''), for any 
A > Ar^^ Also, for any A < A^, (1 + gii)|hgw2(A)p < |hf2W2(A)|2. It is because |hf2W2(A)|2 is 
concave in A and attains maximum at A2^^^ whereas |h^W2(A)p is linearly increasing with A. Since 
the intersection point Ag < A™*"*, we have (1 + gii)|h|^W2(A)p < |hf^W2(A)p which violates the 
constraint set requirement. 

2) Eliminating non-sum-rate optinfial solutions: Now, we combine the above results. For any 
solutions bj G Q.^, we have uj G ^^^ and u £ VL^^ for 6 = 0, 1. This is because uj is on the boundary 
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separating ^^^, ^^^. Thus, if uj is sum rate optimal in ^^^ but not sum rate optimal in ^7^^, then the 
sum rate optimal solutions in ^7^^ achieves a higher sum rate than uj which is in the same constraint 
set. Then, u can be removed from the candidate sets of the maximum sum rate point over constraint 
sets Q^''[J^^^. Applying this argument, we combine the following: 

For any u € Q^ which maximizes Z4 in ^^^ is also in il^*' and achieves a smaller Z3 than other 
sum rate optimal solutions in 17^'^, namely W2(A2 ) U |TTr\- Thus, we have: 

Ifcj = arg max min (Zg, ^4} , then w G wafA^) I J — ^ I J V^'''^ (95) 

waenioUs^oo '^ J ^ ||n22|| ^ 

Similarly, for any uj £ il^ which maximizes Z2 in the constraint set Q,^^ is also in fl^^ and achieves 
a smaller Z4 than other solutions. Thus, we have: 

If w = arg max min (^2,^3,^4}, then w G W2 (A^) I J — ^ I J Vo '*. (96) 

waenoiUnioUn™ '^ J ^ ||n22|| ^ 

The candidate set remains unchanged because Wj performs worse than other solutions and itj^^ G 

Vf. 

Lastly, for any ujh £ U^ which maximizes Zi in constraint set Q^^ is also in Q^^ and achieves a 

smaller Z3 than other solutions. And for oJq G il which maximizes Zi in constraint set Q^^ is also 

in Q^^ and achieves a less Z2 than other solutions. Thus, we have: 

If (J = arg max min f Zi, Z2, Z3, ^4} , then w G W2(A^) I J — ^ I J 1)2'''^. (97) 

W2enii|Jn"i|Jnio|Jnoa '- J ^^ I|h22|| 

Therefore, the final candidate sets are 

This result holds for wi because the optimization problem is symmetric. 

To compute the closed form A2 where W2(A2 ) balances channel powers, we can use the same 
approach as before to obtain A^ = ttt — ,12 , /-, — ■:,-u'^'^u2 L-uh^ — ,-, , 

F. Proof of MRT optimality conditions in the ND region 
Define the following two functions in < A2 < 1: 

• i*'i(A2) is a concave function in A2 and attains maximum at A^^'^. 

• ^2(A2) is an arbitrary function in A2 but satisfy the following properties: 

- for any A > A^™^ F2(A) < F2(A^«^). 
<0. 



dF2{X) 



d\2 



A2=A^»RT 
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Lemma 8: The condition i^2(A2™^) > Fi{X2^'^) is a necessary and sufficient condition for 

A^^^ = argmaxmin{Fi(A2), F2(A2)}. (99) 

0<A2<1 

Proof: 
"=>": Let A2 = argmaxQ<;^2<;^ min{Fi(A2), -F2(A2)}- Denote a set A 

A = {A2:Fi(A2)<F2(A2)} (100) 

Let A2 G A and therefore i^i(A2) < ^i(A2) < Fi(A2™^)- The second inequahty is due to the 
fact that Fi(A2) attains maximum at A^'^'". Note that i^i(A^^^) < F2(A^^'') by assumption, 
thus A^^"^ G A and we can write 

Fi(Ar^) < Fi(A*) < Fi(A^'^^). (101) 

Since -F'i(A2) is a concave and has unique maximum at A2''^'^, we have A2 = A2''^'^. 
"<^": We start with A2™^ = argmax;^ min{Fi(A2),-F2(A2)} and we proceed with contradiction. 



9^2 (A) 



Assume Fi(A^'^^) > F2(Ar^). Since -^ 



< 0, there exist A, = A,™'^ - e with 

A2=A^^RT ' ^ ^ 



arbitrary sniaU e > 0, which satisfies -Fi(A2 ) > -F2(A2 ) > -F2(A^^'^) and contradicts to the 
assumption that A2"^'^ = argmax;^^ min{Fi(A2), -F2(A2)} • 

■ 
Note that the sum rate in the ND region is 

R"'^= max min{C2(A2)+Ti(Ai,A2),C2(A2) + A(Ai,A2)}. (102) 

0< Ai,A2<l 

From Thm. [21 we can write the Pareto optimal beamforming vectors w^ in the ND region as a 
function of the real valued parameter A^ in (|21|) . Hence, we can rewrite the rate expressions in ^ 
as functions of A,-, i = 1, 2: 



2C2(A2)+T,(A„A2) = l+52(A2)P™.. + Ai||h2ifP, 
2C2(A2)+i..(A,A2) ^ (1 +52(A2)P...) fl + — ^ 

V 1 + A2 h 



gi{Xl)Pma. ^ (^°^) 



12 1 1 ^max 

where gi(Ai) = (\Ai||n2ihii|| + ^1 - Ai||n^ihii||) and5r2(A2) = (\A2lln12h22ll + VI - A2||n^2h22||) ■ 
Since logarithm function is monotonic, it does not change the maximization solution and from now 
on, we consider maximizing the minimum of the following two functions, 

Fi(Ai,A2) = l+52(A2)P™a. + Ai||h2lfP„,, (104) 

ff i(Ai)P„ 

*^2||hl2||"J'maa;. 



F2(Al,A2) = (l+ff2(A2)P^ax) l+ ./;\;^ "p - (105) 

I + A2 hi2r>V 
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Lemma 9: It can be shown that the function gi{Xi) is concave in A^ for i = 1,2. Fi(Ai,A2) is 



concave in Ao and attains its maximum at X^^"^ = .J ?i2iiu 112 - 

^ ||h2i|K||hii||^ 

Proof: Note that the first and second derivatives of gi{\i) with respect to A^ are -^gi{X 



in--h'-iP - iin-^h'-iP + iin-'h'-iiiin-^--h-'ii (—^ 



2Xi 



and #75i(Ai) 






dXi 

Since A,' is 



between zero and one, the second derivative of gi{Xi) is always negative, for all A^. Set the first 
derivative to zero and we obtain the maximum A^^'^ = .J 112 nu 112 - ■ 

^ ||h2i||-^||hii||^ 

Lemma 10: F2(Ai,A2) satisfies 

9i^2(Ai,A2) 



dXo 



<0. 



Proof: 

^/2{X,A2) 



\2=\l 



-^ [l + g2{X2)P„ 



+ {l+g2{X2)Praax)(-^{l + 



gi{Xi)P„ 



l + A2||hi2PP™ 

gi{Xl)Pmax 
1 + A2||hi2|pPmax 



X ^;^MRT 



(a) 



^1 + g2{X2)Pmax) 



(JL(, 



gi{Xi)P„ 



1 



i22 



^P 

1 rr 



max 

< for any < Ai < 1 



Vc>A2 V l + A2||hi2PP„ 

gl(Al)||hi2|pPmaa; 

1 _l_ \MRT||v, 112 p ^2 
i + Ao ||ni2|| ^max) 



I MRT 



where (a) is due to the fact that g{X2) is concave and attains its maximum at A2 

Lemma 11: For any A > A^™^, F2(A) < F2(Ar"). 

Proof: For any A, 52(A) < 52(A2™^) because g2{) is a concave function and attains maximum 
at A^^'^. Also, for any A > X^^'^ , the denominator of -F2(A2) increases. Thus, for any A > Ag^^'^, 
i^2(A) < F2(Ar"). ■ 

By Lemma [HI R"''^ is maximized by A^'^'^ for arbitrary fixed Ai if and only if Fi[Xi,X2^'^) < 
i<2(Ai, A2^^'^) which is equivalent to the following: 



i^i(Ai,Ar"^)<F2(Ai,A^ 



MRT ^ 

2 ) 



^ Alllh; 



21 



< 



l+||h22fP,. 



l + ||hi2||2cos2(02)^„ 



-ffi(Ai). 



(106) 



Also, for arbitrary fixed A2, Fi(Ai,A2) is linearly increasing with Ai and i^2(Ai,A2) is concave in 
Ai and attains maximum at Af^'^. Similar to the argument before, there are at most 2 intersection 
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points between Fi(Ai, A2) and F2(Ai, A2). We observe that 

< (1 + 92iX2)P^a.) (1 + TT¥§rWF^) (10^) 

\ i + A2||ni2|| r^max/ 

= i^2(0,A2). 

Note that gi{0) > except when hn is orthogonal to h2i whose probability is zero almost surely. 
Since Fi{0, A2) < -^2(0, A2) for any A2, there is at most 1 intersection point. If there is no intersection 
point, the curve F2(Ai, A2) is above Fi(Ai, A2) for any Ai and therefore the optimal value of Ai which 
maximizes Fi(Ai, A2) at Ai = 1. If there is 1 intersection point, denote the intersection solution as 
A^ . Graphically, it is clear to see that if and only if A^ < A^^'^, the optimal solution is Ai = A^^^'^. 
Thus, we have for arbitrary fixed A2, 

Ai = 1 is optimal if and only if -^1(1, A2) < i^2(l, A2) 
Ai = Ar^ is optimal if and only if A^ ^ < A^^ 

where A^ is given in (|71|) . Now we combine the conditions for Ai and A2 and after some manipula- 
tions, we obtain the following: 

(A^^^, A^^^) is optimal if and only if 

Cllln^lhlllP ,,„,2.. ^^ (l + ||h22fP^a.)||hn|P ,.^^. 



C2||h2iP-2^/3I^|h^ihn|+ci||hnP ' ' ' (1 + ||hi2P cos2(^2)P™a.)||h2i 

(Ai = 1, Ar^) is optimal if and only if ||h2if < (1 + ||h22f P^a:.)^ irJ^T^^^^^TT^TT^^) 



l + ||hi2||2cOs2(02)P„ 

hl2|pCOS^(82)Pmt.x + l 



where Ci = uy,,,r2^J!T(Z\P — XT and C2 - ^' 



G. Proof of MRT optimality in the DD region 

In this section, we provide the sum rate optimality conditions for two MRT strategies, namely: 
interference amplifying beamforming w^ = ,,^^',, , in Section IXI-G.ll and direct channel beamforming 
Wj = 11^, in Section IxTGH 

1 1 -til i II 

1) Optimality conditions of amplifying interference in the DD region: We aim to prove that the 
beamforming vector w^ = uy^'\, is sum rate optimal in R'^'^ if and only if (1 + (7jj)||hii|p cos^(0i) > 
\\hjif and \\hii\\^ cos^ [6 i)P„,^^ < || - 1, where gkm = ||h|^„w„j|p. 

Due to symmetry of the problem, the proof for wi and W2 is similar and we only give the proof for 
W2 here. First, by the definition of Q'^^ (IMl) . we observe that W2(A2*^'^) = t^^^ is in the constraint 
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set Q^^ if and only if the following constraints are satisfied: 



vl+5ii 



hgw2{XTl >|hf2W2(A 



MRT \ I 2 



p <^-l 

-* max ^ ^ 



9ii 



which are equivalent to 



hgw2(Ar") 

l + 5ll)l|h22f COS^(^2)>||hl2f 

a ■ (110) 

\\h22fcOs\e2)P^a.<—-l 

911 

Thus, if and only if (|11(J|) is satisfied, the beamforming vector W2(A2™'^) = ijif^ is in il°^. Now 
we establish that this is the sum rate optimal solution. 

Zi is monotonically increasing with g22 in constraint set Q^^. Note that W2(l) = ttjtti D is not in 
the constraint set Q^^. We see this by observing the constraint set of Q^^ requires: 
. |hgw2(A2)P = A2||h22|P . Thus, we have |hgw2(Ar^)|2 < |h^2W2(l)p. 

• |hf2W2(A2)P = (\/A2||n22hi2|| + \/l — A2 1| 11^2^12 1|) IS coucavc in A2 and attains maximum at 
A^^ where W2(Ar'') = p^- Thus, |hf2W2(l)|2 < |hf2W2(Ar'')P- 
If the conditions in (IllOp are satisfied, we have 

(l+5ll)|hiw2(l)|2 > (l + <7ll)||h22f COS2(02) > ||hl2f > |hf2W2(l)p. 

To satisfy the constraints of both fi"^ and fi^^ the sum rate optimal solution lies on the boundary 
between Q^^ and ^^^, namely i}^. 

The constraints set Q,^^ is empty. Using the same argument as in the case of Zi, for any A2 that 
satisfies |h^2W2(A2)P > ^ - 1 must satisfy A2 > A^'^^. Also, any A2 > A^'^^ satisfies 

(l + 5ii)|hiw2(A2)P > (l + 5ii)l|h22fcos2(e2) > ||hl2f > |hf2W2(A2)p. (Ill) 

Thus, for any A2 that satisfies B = must have A = 0, whih indicates that the constraint is empty. 

Similar to the argument before, |h|^W2(A2)P > ^^ — 1 is a tighter constraint for |h22W2(A2)P than 

(l+5ii)|h^W2(A2)p > |hf^W2(A2)p in 17'^°. As Z4 is monotonically decreasing with |h^W2(A2)p and 

increasing with |hf^W2(A2)P in ^^^ , the sum rate optimal solution in this case is the beamforming 

vector which satisfies: 

|hiw2(A2)P = ^-l 

511 

:i+5ll)|hf2W2(A2)P = |hf2W2(A2)P 
^From now on, we write W2(A2 — 1) as W2(l). We must not confuse this with the first element of vector W2. 
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which is in Q^ . 

Now we show that W2(A2'^'^) = uy^ is the optimal solution in Q^^. Z2 is monotonically increasing 
with 512. Since we assumed that t^^^ is in Q'^^, it is the optimal solution. 

Finally, we notice that for any sum rate optimal solutions in ^^ which maximizes Zi in ^^^, it 
is also in the constraint set Q^^ and therefore achieves a smaller sum rate than ttt-^. Similarly, any 
solution in il^ that maximizes Z4 is also in constraint set 57°^ and therefore achieves a smaller sum 
rate than Trrr^- 

||hi2|| 

2) Optimality conditions of direct channel beamforming in the DD region: Now we prove that the 
beamforming vector Wi(Ai) = m. "m attains the maximum sum rate in the DD region, for arbitrary 
fixed w, if (1 + 5,,)|hfw,(l)|2 < |hjfw,(l)|2. 

We provide the proof for W2 for simplicity as by reversing the role Tx 1 and Tx 2, the proof for wi 
can be obtained. Notice that if the optimality condition is true: (1 +5'ii)|h22'W2(l)P < |hi2'W2(l)P, 
then the following arguments are true. 

The constraint sets i7°^ and 17°" are empty. This is because 



h: 



H n22hi2 



"ilonJll 



2 



< |hf2W2(0)| 



2 



l + 5ii)|hiw2(0)p = (1 + 511 

22"12 

Thus, together with the assumption above: (1 + 5ii)|h^W2(l)P < |hf^W2(l)P, we have 

f |hf2W2(l)P > (1 + gn)|hiw2(l)|2 

[ |hf2W2(0)p > (l + gii)|hf2W2(0)|2. 
Since |h^W2(A)p is linearly increasing with A and |hf^W2(A)p is concave in A, we draw the conclusion 
that for ah A, |hf2W2(A)p > (1 + 5ii)|hf2W2(A)|2. Thus, for ah \, A=l. 

The sum rate optimal solution in ^^-^ is either Q^ or Tr^f^- If the optimal solution is tt^^ then 
we know that Q^^ is empty and thJ^ is sum rate optimal. If the sum rate optimal solution in Q^^ 
is in fl^, then these solutions are also in constraint set Q^^ which achieve a smaller sum rate of Z^ 
than j^. 

||h22|| 

Now, we obtained the MRT optimality conditions for each transmit beamformer w^ given Wj. 
Apply the same approach and reverse the role of Tx 1 and 2, we obtain the conditions for Wj. 
Combine both inequalities to obtain the conditions as shown in Theorem [71 

References 

[1] Z. K. M. Ho, D. Gesbert, E. Jorswieck, and R. Mochaourab, "Beamforming on the MISO interference channel 
with multi-user decoding capability," in Proceedings of Asilomar Conference, 2010, pp. 1-6. 

46 



[2] A. B. Carleial, "A case where interference channel does not reduce capacity," IEEE Transaction on Information 

Theory, vol. IT-21, no. 5, pp. 569-570, 1975. 
[3] M. H. M. Costa, "On the Gaussian Interference Channel," IEEE Transaction on Information Theory, vol. IT-31, 

no. 5, pp. 607-615, 1985. 
[4] X. Shang, B. Chen, and H. V. Poor, "Multi-User MISO Interference Channels with Single-User Detection : 
Optimality of Beamforming and the Achievable Rate Region," submitted to IEEE Transaction on Information 
Theory, 2009. 
[5] A. S. Motahari and A. K. Khandani, "Capacity Bounds for the Gaussian Interference Channel," IEEE 

Transaction on Information Theory, vol. 55, no. 2, pp. 620-643, 2009. 
[6] I. Sason, "On Achievable Rate Regions for the Gaussian Interference Channel," IEEE Transaction on Information 

Theory, vol. 50, no. 6, pp. 1345-1356, 2004. 
[7] R. H. Etkin, D. N. C. Tse, and H. Wang, "Gaussian Interference Channel Capacity to Within One Bit," IEEE 

Transaction on Information Theory, vol. 54, no. 12, pp. 5534-5562, 2008. 
[8] Y. Weng and D. Tuninetti, "On Gaussian Interference Channels with mixed interference," in Proceedings of the 

2008 Information Theory and Applications Workshop (ITA 2008), 2008, pp. 1-5. 
[9] T. S. Han and K. Kobayashi, "A New Achievable Rate Region for the Interference Channel," IEEE Transaction 
on Information Theory, vol. IT-27, no. 1, pp. 49-60, 1981. 
[10] H. Sato, "The capacity of the Gaussian Interference channel under strong interference," IEEE Transaction on 

Information Theory, vol. IT-27, no. 6, pp. 786-788, 1981. 
[11] V. R. Cadambe, S. A. Jafar, and S. Shamai, "Interference Alignment on the Deterministic Channel and 
Application to Fully Connected Gaussian Interference Networks," IEEE Transaction on Information Theory, 
vol. 55, no. 1, pp. 269-274, 2009. 
[12] S. A. Jafar and S. Vishwanath, "Generalized Degrees of Freedom of the Symmetric Gaussian K User Interference 

Channel," IEEE Transactions on Infromation Theory, vol. 56, no. 7, pp. 3297 - 3303, 2010. 
[13] B. Bandemer, A. Sezgin, and A. Paulraj, "On the Noisy Interference Regime of the MISO Gaussian Interference 

Channel," in Proceedings of Asilomar CSSC 2008, 2008. 
[14] S. Vishwanath and S. A. Jafar, "On the Capacity of Vector Gaussian Interference Channels," in Proceedings of 

Information Theory Workshop, 2004. 
[15] M. Charafeddine, A. Sezgin, and A. Paulraj, "Rate Region Frontiers for n-user Interference Channel with 

Interference as Noise," in Proceedings of Allerton Conference, 2007. 
[16] E. A. Jorswieck, E. G. Larsson, and D. Danev, "Complete Characterization of the Pareto Boundary for the 

MISO Interference Channel," October, vol. 56, no. 10, pp. 5292-5296, 2008. 
[17] J. Lindblom, E. G. Larsson, and E. A. Jorswieck, "Parameterization of the MISO IPC Rate Region : The Case of 
Partial Channel State Information," IEEE Transactions on Wireless Communications, vol. 9, no. 2, pp. 500-504, 
2010. 
[18] S. G. Kiani and D. Gesbert, "Distributed Power Allocation for Interfering Wireless Links Based on Channel 
Information Partitioning," IEEE Transactions on Wireless Communications, vol. 8, no. 6, pp. 3004-3015, 2009. 
[19] A. B. Carleial, "Interference Channels," IEEE Transaction on Information Theory, vol. IT-24, no. 1, pp. 60-70, 
1978. 



47 



[20] R. Mochaourab and E. A. Jorswieck, "Optimal Beamforming in Interference Networks with Perfect Local 

Channel Information," will appear in IEEE Transaction on Signal Processing, pp. 1-30, 2011. 
[21] Y.F. Liu, Y.H. Dai, and Z.Q. Luo, "Coordinated beamforming for MISO interference channel: Complexity 

Analysis and efficient algorithms," submitted to IEEE Transaction on Signal Processing, 2010. 
[22] V. S. Annapureddy and V. V. Veeravalli, "Sum Capacity of MIMO Interference Channels in the Low Interference 

Regime," submitted to IEEE Transactions on Information Theory, 2010. 
[23] D. Tomecki and S. Stanczak, "On Feasible SNR region for multicast downlink channel: Two user case," in 

Proceedings of IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP'IO), 2010, 

pp. 3474-3477. 



48 



